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Abstract 

We modify the time-dependent electric potential of the Paul trap from a 

sinusoidal waveform to a square waveform. The exact quantum motion and 

the Berry's phase of an electron in the modified Paul trap are found in an 

analytically closed form. We consider a scheme to detect the Berry's phase 

by a Bohm-Aharonov-type interference experiment and point out a critical 

property which renders it practicable. 
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I. INTRODUCTION 



The Paul trap is an instrument which suspends free charged and neutral particles without 
material walls. Such traps permit the observation of isolated particles, even a single one, over 
a long period of time The Hamiltonian of the Paul trap has the form of a time- dependent 
harmonic oscillator, 

Hit) = ^p\t) + lmu;'{t)q'{t) (1.1) 

whose effective spring constant is of the form [Q] 

k{t) = a + bcos{2'Kt/T). (1.2) 

The quantum motion of the Paul trap has been studied in Refs. . It is well known that 
the generalized invariant for Eq. ( p..lD can be written as 

m = 9-{t)j + 9o{t)^^^ + 9At)j. (1.3) 
Here, using the classical solutions satisfying 

A,2(t)+^'(t)/l,2(t) = 0, (1.4) 

we have 

g^t) = Cflit) + C2h{t)h{t) + C3/|(t), 

(?o it) = -m{c^j\{t)f\{t) 

+ {c2/2)[h{t)m + mMt)] (1.5) 

+C3/2(t)/2(t)}, 

g+{t) = m'[cJ!it) + C2h{t)m + C3/|(t)] 

where Ci, C2, and C3 are arbitrary constants. 

Recently, Ji et al. found the exact eigenfunctions of I{t): 
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where if„ is Hermite polynomial. For a time-periodic quantum harmonic oscillator, analyz- 
ing the wave function in Eq. (|1.6| ), they constructed a cyclic initial state (CIS) such that 
^^(t + r') = e^'^"(^'Vn(^) with 

Xn(r') = -(^ + ^) r' ^^,dt. (1.7) 
V 2/ JO mg_[t) 

and calculated the corresponding Berry's phase (see Ref. for the Berry's phase and Ref. |§] 
for its nonadiabatic generalization). Subsequently, a new type of CIS, whose period is a 
multiple of the period of the Hamiltonian, was found 0. 

In this paper, we modify the time-periodic electric potential from the sinusoidal waveform 
in Eq. ( |1.2|) to a square waveform. This square potential has stable classical solutions, as 
the sinusoidal potential does. This means that we can suspend charged particles using this 
modified potential, as we do in the original Paul trap. Furthermore, the classical solutions of 
this modified Paul trap are very simple, so we can calculate the exact quantum solutions in 
a simple closed form. (Note that the classical solutions of the original Paul trap are Mathieu 
functions, which are difficult to deal with.) The purpose of this paper is to find the Berry's 
phase for the modified Paul trap and to propose an experimental scheme to detect it. 

As seen from Eq. ( |1 . 5| ) , there is an arbitrariness in fixing the invariant, and hence the 
complete set of the Fock space (eigenstates of the invariant). Therefore, we should show 
that the phase change of an eigenstate, Eq. ( p.. 61) , is irrelevant to which invariant we choose. 
There is another problem: When we let the electron beam pass through the modified Paul 
trap, it seems that we should have a single eigenstate for a coherent interference pattern. 
However, it turns out that if we prepare a plane wave of the electron - which can be expanded 
as the eigenstates in Eq. (|1.6| ) - we get a coherent interference pattern. 

In Sec. we apply the result of Refs. 0], 0, and [^] to the modified Paul trap to find the 



exact quantum state and the Berry's phase. In Sec. |T|, we present a Bohm-Aharonov-type 



experimental method for detecting the Berry's phase of this system. The key feature which 
renders this experiment practicable is that the phase change is independent of the invariant 
we choose, and for a coherent interference pattern it is sufficient to prepare a plane wave 
entering the trap. A summary and discussions are given in the last section. 



II. EXACT QUANTUM MOTION OF THE MODIFIED PAUL TRAP 

A. Quantum Mechanics of the Paul Trap 

The classical and quantum motion of an electron in the Paul trap is described by the 
following Hamiltonian [jl|: 

H{t)=H,{t) + Hy{t) + H,{t) (2.1) 

where 

Hx = + ^mulx'^, (2.2a) 

Hy = + Imuy, (2.2b) 

= (2.2c) 

Here, the Hamiltonians of the x- and the ?/-motions have the form of a time- dependent 
harmonic oscillator with 

where 

$(t) = U + V cos(27rt/r). (2.4) 

is an applied voltage, d is the gap of the walls of the Paul trap, and e is the absolute value 
of the electron charge. 

The wave function of this system satisfies the time- dependent Schrodinger equation 
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y, z, t) = H{t)^{x, y, z, t). (2.5) 

Using the method of separation of variables, we have three independent equations: 

t) = Hi{t)^i{ii, t), (z = X, y, z). (2.6) 

Here, the equation in the z-direction gives the plane- wave solution \E'^(z,t) = Q'^i^zz-E^t) ^ 
addition, since Eqs. ( |2.2a| ) and ( p.2b| ) are the Hamiltonian of a time-dependent harmonic 
oscillator, we can find and using the methods found in Refs. Q and 0. 



B. Modified Paul Trap 

Now, we modify the applied voltage from the form in Eq. ( p. 41) to the following square 
wave form (see Fig. |I]): 

$1 > 0, for T2 — T < t — rr < —T2, 

m = { (2.7) 

$2 < 0, for -T2 <t-rT < T2, 

where r is an integer. Then, the frequencies of and Hy are described by 

Uf, for T2-T <t~rT < -T2, 

ul{t) = { (2.8a) 
-iU2, for — T2 < t — rr < T2, 

, -ujf, for To — r < t — rr < —To, 
col{t) = \ (2.8b) 

UJ2, for — T2 < t — rr < T2, 
where 

2 e|$i| 2 e|'^2| 
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^(t)t 



^1 







2Xi 



2x, 



FIG. 1. The time-dependent potential of a square waveform. 

In order to study the quantum mechanics of this system, it is necessary to determine 
independent classical solutions of Eq. (|L^) for the x- and the y-components with Eqs. (p.8aj) 
and ( p.8b| ), respectively. These solutions are fully analyzed in Ref. 0. Since the effective 
spring constant of the modified Paul trap alternates between positive and negative values, 
we should check carefully that the solutions of Ref. are applicable in this model. After 
tedious calculations, we verified that our classical solutions are identical with the solutions 
of Ref. 1^ with the replacements of uj2 by —iuj2 in the x-component and uoi by —iuoi in the 
y-component. 

As a result, we find the classical solutions for the x-component to be 

^x.rC^^i^*"''^) + S:^.^e-*^i(*-^^), ioi T2-T <t-rT < -Ts, 
C^,re^2(t-rr) _^ ^,e-'^2(*-^^) , for -T2<t-rT < Ta, 

where 



(2.10) 



Ir \ 






2 

V 



y Bx,r J 



(2.11) 



The coefficients Ar^^r and B^ r, belonging to successive values of r, can be related by a matrix 
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P obtained by imposing continuity for fx{t) and its derivative &X,t — —T2+rT and t — T2+rT. 
These lead to 



(2.12) 



with 



V 



(a ] 


pr 


(a \ 




y ^x,r J 










-iPx 



i(^x,2e 



■lUJlT 



JIT ^ 



where 



Oix,i — COS 2ujiT2 cosh 2UJ2T2 H — sin 2ujiT2 sinh 2uj2T2, 
l3x,i — sin 2a;iT2 cosh 2a;2T2 — - cos 2a;iT2 sinh 2uj2T2, 
/3x,2 = |sinh2a;2T2, 



and 

LO\ UJ2 ^1 ^2 

e = h— , r?= , 

UJ2 i^i CO2 ^1 

where a^^i, (3x,i-i and 13x^2 satisfy the condition 

Solving the eigenvalue problem for the matrix P, we find the eigenvalues 



p± = A. ± JXl - 1 



where A^; = a^^i cosujiT + [3x^i sina;ir, and their corresponding eigenvectors 



oc 



I3x,2e"^'' 
v^, ± i./AE - 1 



(2.13) 



(2.14a) 
(2.14b) 

(2.14c) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



where Vx — eta;,! sina;iT — cosa;iT. If |Aa;| < 1, p± are complex conjugates. Investigating 
the form of the matrix P, it is easy to find that the solutions corresponding to two eigenvalues 



are also complex conjugates. Therefore, the two independent solutions are taken to be the 
real and the imaginary parts of one of them. In this case, we can set 



(2.19) 



with one eigenvalue 



P+ 



(2.20) 



where tan^j; = Jl — XI/ Xx- Then, the classical solution for lA^^I < 1 can be written as 



3'"^^ [A^,oe''^i^*~''^) + S^,oe-^'^i(*-''^)] , for T2 - T < t - rr < -T2, 

2''''^" [Cx,Oe^^^^~''^^ + L'^,oe-'^2(t-rr)J ^ ior -T2 <t-rT <T2. 

In the same way, we have the classical solution for the y-component: 



fyit) 



Ay^oe^^('-'^^ + S^,oe-'^i^*-''^)] , for ts - r < i - rr < -T2, 
Cy,oe*'^^(*-'^^) + Dy^oe-'''^^^-'-^^] , for - T2 < t - rr < T2, 



with 



(2.21) 



(2.22) 



(2.23) 



^ Cyfl \ 1 (1 - M/^2)( 



-(jJlT2+iu)2T2 



[1 + iLUi/iU2)e' 



UJlT2+i0J2T2 



^ (1 + iLJi/LJ2)e-'''^^-'''^^^ (1 - iu}i/u}2)e'^'^^-''^^^^ ^ 



(a ^ 

\^y,o J 



where 



(2.24) 



Uy = ay^i sinha;ir — Py^i cosha;ir. 



(2.25) 



Oiy,i — cosh 2a;ir2 cos 2u;2T2 — - sinh 2a;ir2 sin 2ci;2T2, 
Py^i — sinh 2a;iT2 cos 2a;2T2 — - cosh 2a;iT2 sin 2a;2T2, 
f^y,2 = |sin2u;2r2. 



(2.26a) 
(2.26b) 
(2.26c) 
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and 



where \y = ay^i coshtuir — Py^i sinhcjir. 

The two independent real solutions and f2(t) are given by 



(2.27) 



(2.28) 



for the X and the y components, respectively. These solutions exhibit stable motions for 
A < 1 (A stands for Aa; or A^); that is, they oscillate with bounded amplitudes. 
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FIG. 2. The stability- instability diagram. The vertical strips stand for the stable regions in the 
x-motion, the horizontal strips for the y-motion. 

It is important to know the stable regions in uji-uj2 diagram where the classical solutions 
are stable. For ri = T2, we present the stable regions and the unstable regions in Fig. |^. 
This map is similar to the stability diagram obtained from the Mathieu equation. Only the 
overlapping regions of the x-stable and the ^/-stable regions are of our interest. Therein, 
the motion is stable both in the x-direction and the y-direction. On the other hand, when 
|A| > 1, the solutions diverge at t ^ 00 or t — > — 00 as discussed in Ref. 0. 
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Now, we fix the generalized invariant by fixing ci, C2, and C3 in Eq. ( p.. 51) . For example, 
to find I{t) such that I{t) = i/(to) in the x-component (to denotes the initial time), we fix 
those three parameters as [|1^] 

Px2 + Bio- 2/3x,25^,o cosa;ir 



Cx,l 



m 



{Ph - B 



2 

x.O 



2 



-A(3^^2Bx,osmuJiT 
Cx,2 = -. -2-, (2.29) 

^{Ph - Bio) 

^ _ + Bio + 2/3x,2^ x,0 COSO^iT 



In this way, we can get the exact wave function of the modified Paul trap, Eq. ( |1.6| ), and 
the phase change (which includes the Berry's phase) for a period, Eq. ([L7|). 



III. EXPERIMENTAL METHOD OF DETECTING BERRY'S PHASE 

A. A^r-periodic Wave Function 

In this section, we present an experimental method to detect the effect of the Berry's 
phase. The existence of the CIS is provided by the periodic classical solutions. As discussed 
in Refs. and 0, if it holds in Eqs. (^) and (p:22|) that 



(j) = ^27r (/, N' = integers and N' 0) (3.1) 

(where the index of x and y is understood), the classical solution is A^'r-periodic. Then, 
g-{t) is A^'r/e-periodic (e = 1 for odd A^', e = 2 for even A^'), accordingly, so is the wave 
function in Eq. ( |1.6| ). 

When we have two independent real classical solutions, say fi(t) and f2(t), we can always 
construct the complex solution as 

m = d^f,{t) + {d2 + ^d,)f2{t), (3.2) 

where di, d2, and d^ {did^ 7^ 0) are real parameters. This solution can be written in polar 
form []TI| : 
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m = \m\e'"^'^ (3.3) 

where 

If we set ci = dl, C2 = 2did2, and 03 = ^2 + dl, we have g-(t) = |/c(t)p. Then, the quantum 
phase, Eq. (|1.71) , of the n-th eigenstate, which is also a CIS with a period r', can be rewritten 



as 



Xn{r') = -{n + l-]e{T'). (3.5) 



2 

Now we are ready to prove that the quantum phase in Eq. (p.5|) is independent of the 
choice of the invariant. That is, the phase change of the eigenfunction of the invariant does 
not depend on what values of q we choose. The proof is as follows: If we assume that the 
phase change of the classical solution in Eq. ( |3.3| ) is altered by varying the parameter values 
di{i = 1,2,3) or q, then there are many classical solutions corresponding to the respective 
periods. However, this contradicts the fact that the classical solution of Eq. (|1.4D has only 



two independent solutions and that they are the real and the imaginary parts of Eq. ( p. 31) . 
This completes our proof. 



B. Experimental Setting 

Now let us consider the experimental arrangement. Suppose we have a single coherent 
electron beam which is split into two parts, and suppose each part is allowed to enter the 
modified Paul trap, as shown in Fig. ^ 
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FIG. 3. The schematic diagram of the experimental arrangement. 



This experiment is similar to the experiment illustrated by Aharonov and Bohm in 



Ref. |jT2[. After the beams pass through the modified Paul traps, they are combined to 
interfere coherently at the point F. Let us denote the paths A-B-C-F and A-D-E-F by 
path 1 and 2, respectively. The electric potential vanishes in region I so that the wave 
function of the electron is described by a plane wave which propagates along the z-direction: 
^(I) = L-3/2e^^-"-^-* where L is a suitable normalization factor. In region II, the potential 
varies as a function of time according to Eq. ( |2.7| ), and — in the x- and the y- 
directions, respectively, but $i and/or $2 have different values for paths 1 and 2. When the 
electron is in region III, the potential vanishes again. 
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Let and \I/2 be the wave functions that pass through path 1 and path 2, respectively. In 
region I, = = \&2(I)- Then, as they enter region II, the two wave functions suffer 

different potentials in the two different Paul traps. Finally, we have \E'(III) = \E'i(III)+\l'2(III) 
in region III. 

In order to form a sharp interference pattern, it is necessary to have the wave functions 
at point F be of the form 

*(III) = ^(I)e-*®^ + ^(I)e-^®^ (3.6) 

such that the pattern depends upon the phase difference 62 — 61. We emphasize that the 
critical factor which renders our experimental scheme practical is that any plane wave which 
splits at the point A does interfere at F, as required, in the form of Eq. (|3.6|) . 

Firstly, when "^i^ is a single eigenstate of the LR invariant, the phase change can be 
easily obtained by using Eq. (|1 . 7]) , and it is evident the final wave function at F is of the 
form of Eq. ( p.6| ). Next, for a plane wave propagating along the 2;-direction with the wave 
number k, 

w= — = —z (3.7) 
m I 



gives a final wave function of the form of Eq. (|3.6|) , as we will see below. Here, D is the 
length of the Paul trap, and T is a multiple of the minimal period of the CIS, which can 
be controlled by the applied voltage. In this situation, the wave function of the electron 



entering the Paul trap is expanded in terms of eigenstates of the LR invariant [13 



^i..(t)= E C'i,>i,„(a;,t) (3.8) 

n=0,2,4,... 

{"^ly, '^2x, and \l/2y can be expanded in a similar manner.) When the electrons leave the 

Paul trap, using Eq. ( p7?| ) or (p.5|), we have 

00 

<f,,it + T) = E Ci,2n^i,2n(x, t)e-(2"+i)^- . (3.9) 

n=0 

Further, in the phase of each eigenstate, the periodicity of the classical solution means that 
2n6ix{T) is a multiple of 2tt. Therefore, we can write 
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^Ut + T) = e-^^l^(^) Ci,2„^l,2n(x, t) (3.10) 



n=Q 



= e-i^-(^)v[;,,(t). (3.11) 

In the same way, we have 

^i,(t + T) = e-^'i-(^)^i,(t). (3.12) 
Then, we have the total wave function which travels path 1: 

^i(t + T) = e-^[^-(^)+^^-(^)le-^^^^i(t), (3.13) 

for path 2, we have 

Therefore, the phase difference between two paths is 

02 - Oi = ^{[^2.(T) + ^2,(T)] - [eUT) + e^y{T)]} (3.15) 
= ^2(T)-ei(T). (3.16) 

Here, we have omitted the indices x and y since the phase changes for a period are equal in 
the X- and the ^/-directions. 

C. Expected Results 

In this section, we present a typical experimental scheme. In region III, we have the 
detector F, and we have a destructive interference when 

\92-9i\=7i. (3.17) 

This destructive interference of the two wave functions, via path 1 and 2, can be obtained 
by controlling the applied voltage or the velocity of the electron beams. By noting the fact 
that when Eq. ( |3.1|) holds, the phase change over the minimal period is 
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e{T' = N'r/e) = In, (3.18) 

we have two methods to obtain destructive interference. Firstly, we can control the applied 
voltages $1 and $2 so that / = 1 in Eq. ( pTip and N[ = 2A^2 ('where and N2 are the 
values of A^' for path 1 and path 2, respectively). Further, we can control the velocity of the 
electron beam so that 

^(T = N[t). (3.19) 



Then we have 61 = n and 62 = 2tt. Secondly, we can control the voltage values so that / = 1 
and I = 2 with the same N'. From Eq. ( p.l8|) , it is clear that two wave functions interfere 
destructively. In Table 1, we present the numerical values of ui and uj2 for N' = 4,8 with 
/ = 1 and for A^' = 3 with / = 1, 2. 
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TABLES 



Table 1. Numerical values of Ui and uj2 for A^r-periodic CISs (A^ = 2,3,4) 



Fig. I I N' uJiT = UJ2T e{T') 

a 1 4 3.14159 ti/2 

b 1 8 2.30517 7r/2 

c 1 3 2.63690 7r/2 

d 2 3 3.48328 vr 
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The graphs of g-{t) in the x-direction for all the cases in Table 1 are shown in Fig. 4. 
The parameters are fixed as in Eq. ( |2.29| ), and hence the region of g~.{t) = l/m (the figures 
are depicted in m = 1 units) reflects that lit) = H{t), as discussed in Ref. 0. By shifting 
these figures by a half period, r/2, we can also get g-{t) in the y-direction. As expected, 
they are A^r-periodic [N = 2, 3, 4), as are their corresponding wave functions, Eq. The 
original and the shifted figures also reveal that the probability density function |\E'(x, t) p 
spreads in the x-direction and the y-direction alternately. 

^" Fig.4a 

35 r ' 1 ' ' 1 • ' 1 • • 1 ' p^ ' rp 

30 - - 
25 - ■ 




S" Fig.4c 

35 r ' 1 ' ^p ' 1 • ^p ' r 

30 - 
25 ■ 

20 




17 



FIG. 4. The shapes of g^{t) for Table 1. The time is denoted in units of r. 



IV. DISCUSSION 

Modification of the time- dependent electric potential of the Paul trap from a sinusoidal 
waveform to a square waveform gives a quantum solution with a simple mathematical form. 
Therefore, we can verify the existence of the iVr-periodic CIS and propose a method to 
detect the corresponding Berry's phase by experiment. 

We estimate the values of the parameters for a practical experiment. To obtain an 
interference pattern successfully, we should have a Paul trap of submillimeter size {d ~ 
10^'^ m). Considering a Paul trap whose length D is of the order of 10~^ m and an electron 
with a speed on the order of 10^ m/s, r' ~ 10~^ s, we have Ui ~ 10® s~^ from Table 1 and 
$ ~ 1 V from Eq. (|2.9|) . For example, for A^' = 4, / = 1 in Table 1, D ^ 6 cm, d ^ 1 mm, 
and f ^ 5 X 10^ m/s, we have r ^ 6 x 10^^ s, uj(^i^2) ~ 5 x 10® s"^, and |$(i,2)| ~ 1 V. These 
values seem practical for an experimental arrangement. 

There have been many applications and tests of the Berry's phase using an optical 
fiber nuclear magnetic resonance (NMR) [^, etc. |]TB[. Nonetheless, there are no 



experiments about the Berry's phase caused by the quantum motions in phase space. (Note 
that the optical phase effect deals with the phenomena of classical electromagnetism, and the 
NMR experiment investigates the interaction between the spin and the external magnetic 
fields.) Our proposal will be a new experiment to detect the Berry's phase caused by a pure 
dynamics in phase space, and we expect that it will play a significant role in understanding 
the quantum motions of a time-dependent system in phase space. 
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Figure captions 
Fig. 1. The time-dependent potential of a square waveform. 

Fig. 2. The stabihty-instabihty diagram. The vertical strips stand for the stable 
regions in the x-motion, the horizontal strips for the y-motion. 

Fig. 3. The schematic diagram of the experimental arrangement. 

Fig. 4. The shapes of g-{t) for Table 1. The time is denoted in units of r. 
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